Modification of the classical Heisenberg helimagnet by weak uniaxial anisotropy and 

magnetic field 



I. A. Zaliznyak 

P. L. Kapitza Institute for Physical Problems, ul. Kosygina 2, 117334 Moscow, Russia 



CO 

o 

O 

CN 

C 

co 



i 

-4-* 
c3 



M. E. Zhitomirsky 

i. D. Landau Institute for Theoretical Physics, ul. Kosygina 2, 117334 Moscow, Russia and 
Institute for Solid State Physics, University of Tokyo, Tokyo 106 Japan 
(Dated: February 4, 1995) 

Classical ground state of the isotropic Heisenberg spin Hamiltonian on a primitive Bravais lattice 
is known to be a single-Q planar spin spiral. A uniaxial anisotropy and external magnetic field 
distort this structure by generating the higher-order Fourier harmonics at wave-vectors nQ in the 
spatial spin configuration. These features are not captured in the formalism based on the Luttinger- 
Tisza theorem, where the classical ground state energy is minimized under the "weak" condition on 
the length of the spins. We discuss why the correct solution is lost in that approach and present 
an alternative microscopic treatment of the problem. It allows to find the classical ground state for 
general Q for both easy-axis and easy-plane uniaxial second-order anisotropy, and for any orientation 
of the magnetic field, by treating the effect of anisotropy (but not the field) as a perturbation to 
the exchange structure. As a result, the classical ground state energy, the uniform magnetization, 
and the magnetic Bragg peak intensities that are measured in experiment, are calculated. 

PACS numbers: 75.10Hk, 75.25+z, 75.50Ec 
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I. INTRODUCTION 

For more than three decades, helical spin struc- 
tures have been a subject of intensive studies. Having 
been discovered theoretically in the pioneering works of 
Yoshimorii and Villain 2 ", they were found experimentally 
in a large variety of materials. Most extensive work was 
devoted to the investigation and explanation of different 
ordered phases and phase transitions in the rare-earth 
metals from Tb to Tm. These phases were shown to re- 
sult from the intricate distortion of the incommensurate, 
almost ferromagnetic exchange helimagnet by tempera- 
ture and strong crystal field anisotropics £^ However, for 
the reason discussed below, a classical treatment, finally 
generalized by Lyons and Kaplan*, predicts no such dis- 
tortion in the ground state of the quasi-isotropic Heisen- 
berg Hamiltonian at T = 0. 

A simple but nontrivial example of the commensu- 
rate spiral is a triangular magnetic ordering found in the 
hexagonal antiferromagnets of CSN1CI3 type. Their mag- 
netic structure consists of six sublattices, three in each 
hexagonal plane at an angle «120° to each other, whereas 
spins in the adjacent planes are antiparallel. Such com- 
pounds have recently attracted much attention because 
of their pronounced quasi-one-dimcnsional nature (ex- 
change between the adjacent spins along the hexagonal 
axis is much stronger than that between the neighbors 
in the plane) . Experiments^^ have demonstrated large 
deviations of the staggered magnetization and suscepti- 
bility in these materials from the results of the mean 
field (MF) calculations which should be valid for classi- 
cal spins, S 3> 1, and thus the importance of quantum 
fluctuations. Therefore, these are also excellent examples 
of quantum helimagnets which also have been recently 



studied in the literature Aifliii Chubukov 1 - and Tanaka 
et. al. 13 calculated classical spin reorientation process 
in magnetic field (see also Ref. and the antiferro- 
magnetic resonance spectra in the hexagonal antiferro- 
magnets of CsNiCl3 type on the basis of six-sublattice 
model. It appears that even for six sublattices the prob- 
lem is very difficult to handle in this way, and essen- 
tially no explicit analytical results can be obtained for 
the spin-wave spectrum. Our motivation in the present 
work was to develop a simpler and more general proce- 
dure for finding the classical ground-state spin configura- 
tions, which would allow to perform reasonably unsophis- 
ticated spin- wave calculations. We adopt the approach 
similar to that used by Nagamiya 3 and Kaplan 4 . Es- 
sentially, it is extension of the iterative procedure first 
suggested by Cooper et al~ for the easy-plane, almost 
ferromagnetic spirals, to a general helimagnet. How- 
ever, keeping in mind its application to the commensu- 
rate antiferromagnets, we will neither discuss here the 
incommensurability effects, nor changes in the ordering 
wave- vector Q by the anisotropy and magnetic field that 
occur for general Q position. These problems, as well 
as the incommensurate-commensurate transitions, were 
treated elsewhere in the framework of the phenomeno- 
logical approachpiSiii 

We start from the following spin Hamiltonian, which 
describes the low-temperature magnetic properties of a 
large class of crystals where spins of the magnetic ions 
are localized at lattice sites, 

H = J2 Jij SiSj + D ~ 7H ]T S, . (1) 

i,j i i 

Here the first term is the Heisenberg exchange interac- 
tion, the second term describes the lowest-order uniaxial 
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anisotropy, and the last term is the Zeeman energy (7 is 
the gyromagnetic ratio). The anisotropy arises from the 
magnetic spin-spin and spin-orbit interactions. Hence, 
the constant D incorporates a small relativistic ratio 
(v/c) 2 , where v is the average velocity of the electron 
in the atom, and c is the velocity of light. This constant 
is usually considered to be small in comparison with the 
exchange coupling constants Jy which are of the electro- 
static origin. This condition, however, can be violated 
in the quasi-low-dimensional magnetic materials, where 
the superexchange path along some directions is indirect 
and the corresponding constants Jij are very small. Al- 
though such cases, as well as singlet ground state systems 
in which D is intrinsically very large may also be of inter- 
est, in the present paper we will always imply condition 
\D\ <C I I / fulfilled. We also restrict our attention 
to the cases in which summation in Q is performed over 
a single Bravais lattice with N sites. 

In the spin-wave theory, the ground state and the ex- 
citation spectrum of quantum Hamiltonian are ob- 
tained in a semiclassical approximation, based on l/S 
expansion. The starting point for such calculation is 
to find an equilibrium spin configuration with the mini- 
mum ground state (GS) energy Eqs m the classical limit, 
S — > 00. In the case of the isotropic Heisenberg Hamil- 
tonian, a rather complete spin-wave theory describing 
its low-temperature properties has been developedAASiii 
However, the problem becomes extremely complicated if 
the anisotropy which exists in real compounds is included 
together with the magnetic field. In this case, even for 
classical helimagnets, very little general results were ac- 
tually obtained. Some particular cases were successfully 
treated in Refs. I3I and lTH but resulted in rather involved 
and hardly usable expressions (in part, because an at- 
tempt to account for the six-fold anisotropy, in addition 
to the easy-plane anisotropy, made problem much more 
complicated). Here we present calculation of the ground 
state spin configurations for the Hamiltonian ||TJ for dif- 
ferent orientations of the magnetic field, and taking into 
account the anisotropy (but not the field) as a first-order 
perturbation of the exchange Hamiltonian. 

II. FAILURE OF THE WEAK CONDITION 

It is known that ground state of a system of equiva- 
lent classical spins on a simple Bravais lattice, in the ex- 
change approximation, is a flat spin spiral described by 
a wave-vector Q (this includes ferromagnetism (Q = 0) 
and antiferromagnetism (Q = K/2) as particular cases). 
A rigorous proof of this result is obtained by solving the 
mathematical problem of finding the absolute minimum 
of the function QJ, which depends on 3N classical vari- 
ables Sf , under N conditions, 

Sfl = S*,Vi, (2) 

imposed on the lengths of the classical spins. This prob- 
lem is solved by introducing N Lagrange multipliers A,. 



Switching to the Fourier representation to use the lattice 
translational symmetry, one obtains the following system 
of equations for spin configuration which minimizes the 
Hamiltonian 

NJ k S k + e z DS% - Yu X k'Sk-k' = \lH8 kfi , 
k' 

X! S k' S k-k' = S2 ^k,0 ■ ( 3 ) 
k' 

Here J k , S k , and \ k are the lattice Fourier transforms 
of the functions Jy, Si, and Aj, respectively; e z is a unit 
vector along z direction, and 6 k is a three-dimensional 
Kronekker symbol. From Eq. |JSJ we easily obtain the 
classical ground-state energy as 

1 E GS = X S 2 -1hS . (4) 

There are two points to note here. First, Eq. © is a com- 
plicated inhomogeneous system of nonlinear equations, 
for which no general solution is easily found. Second, the 
resultant ground-state energy explicitly depends only on 
the values of Ao and Sq (on the latter only if H 7^ 0), 
so all X k for k seem to be irrelevant. This fact 
encourages one to look for the solution for which the La- 
grange multipliers have the form X k = Ao S k Q and the 
system is largely linearized. Evidently, such procedure 
is equivalent to taking into account only one of the N 
conditions given by the second relation in Eq. or 
to replacing Eq. by the so-called "weak" condition, 
J2 t S 2 = NS 2 . This standard approach was formalized 
by Lyons and Kaplan 5 . Then, the solution is easily found 
in the form of a single-Q spiral, 

S k = S 6 k0 + SqS k Q + S*qS k _Q , 

Ao = NJ Q = Nmm{ J k } , (5) 

which is actually a correct result for D > and H along 
the z axis. In this case one can easily choose vectors Sq 
and Sq in (JSJ that satisfy all N "strong" conditions (J2J 
for any D, or for any direction of the field, if D = 0. Thus, 
according to the Luttinger-Tisza theorem, the resultant 
spin configuration is the solution of 1(3} , ie it minimizes 
the energy (JTJ under the "strong" conditions 

Unfortunately, except for the above important but sim- 
ple cases, simple solution obtained in this way appears to 
be physically meaningless. For example, for the negative 
(easy-axis) anisotropy, the H = magnetic structure pre- 
dicted with the help of the "weak" condition is a collinear 
configuration with spins of varying length, parallel to the 
z axis^ This prediction, in particular, contradicts the 
experimental data in CsNiCi3 and related compounds, 
where the easy-axis anisotropy is known to align the spin 
plane parallel to z-axis and to distort slightly the perfect 
120° exchange structure, initially described by a spiral 
with Q = (i,i,l), in reciprocal lattice units (rlu). In 
fact, the "incorrect" collinear phase, corresponding to 
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the longitudinal helimagnet, also appears in CSMCI3, in 
the narrow temperature interval between the two Tjv**. 

The reason for this contradiction and inconsistency of 
the trick with the "weak" condition is the following. It 
is easy to see, that by weakening the conditions on the 
spin length one may obtain additional minima that do 
not satisfy Eq. ©. These artifacts are thrown away 
by verifying all the conditions explicitly. However, some 
extrema of the initial problem © are also lost in this 
way. This can be visualized by considering the energy 
surface Q in the space of spin components, cut through 
by the "physical" manifold |3J). If the real solution, it the 
global minimum, is one of the intersection points, it will 
be inevitably missed by employing the "weak" condition, 
since it is not an extremum of the function Q in the 
extended space. 



III. TRANSFORMATION TO THE LOCAL 
AXES 



sequent decomposition into a series in Bose operators. 
First, we transform the spins from the crystallografic xyz 
frame to the local coordinate axes xiyiZi. In general, this 
is achieved by two rotations: first, by an angle Oi around 
the y axis, and then by an angle fa around the new x\ 
axis, as shown in Fig. 1 . In the resultant local frame the 
spin operators are 



Si = \\T 2 4 ■ ||T M ||.S 4 = HTiH-S* , (6) 
where the transformation matrices are defined by 



|Ti, 



cos 9i — sin Oi 



1 

sin Oi 

1 

COB fa 

sin fa 





cos Oi 


— sin tj> 
cos fa 



(7) 



In order to perform the spin-wave calculation, one de- 
fines the local coordinate axes, XiyiZi, such that classical 

spin at a site i points in direction. The spin operators TL = Si 

are then transformed to these axes and expanded into ij 
a series in Bose operators, eg using Holstein-Primakoff 
transformation. Because the local axis Zi is a classical 
equilibrium direction for spin at site i, the linear in spin 

deviations (ie, in operators af , hi) terms are absent in the Here the matrices 



The transformed Hamiltonian becomes, 



-,(cx) 



Sj+D^i^ 

i 

-jH^mW^Si. (8) 

i 

and TV; , which define 



resultant boson Hamiltonian. Starting from this point, 
one can develop a slightly different procedure to find the 
classical ground state of the Hamiltonian (JTJ and its sub- 



the bilinear in spin components forms coming from the 
anisotropy and the exchange coupling, respectively, are 
given by 



and, 



sin 2 Oi sin fa sin Oi cos Oi 

sin fa sin Oi cos Oi sin fa cos 2 Oi 
— cos fa sin Oi cos Oi — sin fa cos fa cos 2 Oi 



— cos fa sin Oi cos Oi 

— sin fa cos fa cos 2 Oi 

cos 2 fa COS 2 Oi 



(9) 



-,(cx) 



cos Oij 
— sin On sin< 



COS( 

sine 



i COS ( 
! COS C 



sin 0^ sin fa; 

*J + 
'3 - 



: sin ( 
: sin i 



cost 
sin c 



i sin i 
sin ( 



- sin 0^ cos <pj 
ij — cos Oij sin < 

'j 



; COS( 
COS ( 



(10) 



where Oij = Oi — 0j . The angles Oi and fa must be chosen 
so as to cancel the terms S?Si and SfSj which give rise 



to the linear, in af,a,i, contribution to the Hamiltonian 
(JHJ). This leads to the following system of 2N equations: 



2 Jij cos fa sin Oij = D cos faj sin 20 j + ^H x cos 0j — 7-ff z sin 0j , 

i 

2 (cos fa. sin fa — sin fa: cos fa cos 0^ ) 

i 

= D sin 2faj cos 2 0j — jH x sin fa sin 0j + ^H y cos fa — lH z sin fa cos 6j , 



(11) 
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as a zero approximation and then look for the correc- 
tions, expanding them in powers of D. A very similar 
approach was developed by Andreev and Marchenkc 19 
in the framework of the phenomenological Lagrangian 
theory. They found a rather simple and general solution 
of the problem in the leading order of perturbation. 



IV. EASY-AXIS ANISOTROPY (D < 0) 

First, we consider the more complicated case of the 
easy-axis anisotropy, which leads to the most general dis- 
tortion of the exchange structure. Without the loss of 
generality we can set the magnetic field to lie in the yz 
plane. In this case the "right" unperturbed spiral should 
contain the z axis already at H = 0. Therefore, we look 
for the solution of Eq. I|11J1 in the form, 



FIG. 1: Fransformation to the local axes for the distorted 
exchange spiral. 



Q r i 



sin (pQ 



2(J - J Q ) 



(12) 



where 7 = 7/S. Evidently, this system for the angles 0i 
and cf>i cannot be solved for arbitrary D and H. As it was 
mentioned above, only the cases D — and D > 0, H || 2 
can be treated explicitly 1,11 . However, we can start with 
the Heisenberg exchange Hamiltonian in a magnetic field 



In most cases we will consider only the leading (ie, the 
lowest order in D) corrections 59i, 5<fii. Thus, substitut- 
ing l|12|) into (|llfl . it is sufficient to expand the trigono- 
metric functions to the first order in 59i, 8<pi. As a result, 
we obtain the following equations for the deviations 



2_, 2Jij ( cos 4>a cos(Qr i j)58i — sin <fr sh^Qr^)^) 

i 

= 2Jq cos (j) a 89j + Dcos(j)osin(2Qrj) — •yH z (cos(Qrj)56j + sm(Qrj)) , 
2 Jij [(cos 2 0o + shi 2 </>o cos(Qr i j))S(p i + sin (f> Q cos cj> sin(Qr y )<5^] 

i 

= 2(Jo sin 2 0o + Jq cos 2 4>o)5(j)j + Dsin^o cos0o(l + cos(2Qrj)) 

— {^H y sin </> + 7-Hz cos 0o cos(Qrj))S(f>j + jH z sin <j) sm(Qrj)S8j (13) 

I 



which are rather easy to solve. Indeed, upon convolu- 
tion with cos(Qr i j), the variations S9i and Sifii should 
yield the trigonometric functions cos(Qr 3 ), sin(Qr 3 ), 
cos(2Qrj), sin(2<5r 3 ), etc, which must cancel the right- 
hand side of the Eq. I|ll|l . Thus, they should be searched 
in the form of an expansion, S9i = J2n( a " cos(nQr i ) + 
P n sin(nQr i )), where the order of the coefficients a n , (3 n 
is H n and \D\ n ' 2 . 

An important remark is due here. One can also look 
for the solution of Eq. ifTTJl with 0i — ip + Qr { + 59^ 
instead of Eq. (|12fl . It differs by rotation of the arbi- 
trarily chosen "first" spin, and thus of the whole spin 
structure, by an angle tp within the plane of the spiral. 
Evidently, in the leading approximation (to the first or- 



der in D) all such structures have the same energy, as 
it was the case for the pure exchange. This fact reflects 
the remaining continuous degeneracy of the ground state 
with respect to such rotations. Moreover, by perform- 
ing the corresponding expansions, it can be shown in a 
straightforward algebraic way that this degeneracy ex- 
ists for all orders in H and \f\D\ up to n which satisfies 
the condition nQ = 0, if such n exists. This condition 
means that spin structure is a commensurate spiral and 
can be described in terms of n sublattices. In this case 
the continuous degeneracy corresponding to spin rota- 
tions within the spiral plane is lifted in the order H n and 
|£)|ra/2_ ^phg g rounc j state preserves only n-fold degener- 
acy corresponding to the arbitrary choice of the "first" 
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spin, ie, to the rotation of all spins in the exchange struc- 
ture by an angle {2irm)/n in the plane of the spiral. In 
the case of the incommensurate helimagnet, where this 
discrete symmetry is absent, the continuous degeneracy 
of the ground state is not lifted. 

The above conclusions can be established from very 
general arguments. Evidently, a perturbation of the n- 
sublattice exchange spiral by adding I?^ i cos 2 ('0 + #i) 
and -ff ^ i cos(V' + Qi) to the exchange Hamiltonian is 
invariant with respect to the transformation ip — ► if) + 
(2vrm)/n. Thus, the perturbed ground-state energy 
should not contain harmonics below cos(nV>) and sin(n^). 
Hence, explicit dependence on the angle ip m its Taylor 
expansion with respect to the perturbation parameters D 
and H only appears in the n/2-th and n-th order terms, 
respectively. Another way to establish this fact is based 
on exchange symmetry arguments 19 . Exchange structure 
with n sublattices has an n-fold rotation axis perpendicu- 
lar to the spin plane. Thus, the lowest-order perturbation 
of the exchange energy which has appropriate symmetry 
and which lifts such rotations must be proportional to 
cos(rM/>). In the case of the uniaxial anisotropy which al- 
ready contains the two-fold symmetry, it is achieved in 
the n/2-th order of perturbation. Magnetic field, on the 
other hand, can lift the n-fold rotation axis characteris- 
tic of exchange spiral only in the n-th order of pertur- 
bation. With the continuous GS degeneracy removed, 
the energy gap appears in the corresponding (formerly) 
Goldstone mode. This gap is induced by an "effective 
field" ~ \D\ n / 2 , or ~ H n , and, therefore, is of the order 
~ \D\ n / 4 , and ~ iW 2 , respectively. 

This situation was closely studied in the hexagonal 
CsNiCl3-type antiferromagnets, which are our reference 
in the present paper. It was shown, that in the easy-axis 
case the structure with t/j = 0, described by Eqs. I|12fl . is 
stabilized by the anisotropy in the third order of pertur- 
bation, and that the magnon gap ~ |Z?| 3 / 2 appearsi22i2i 
For the easy-plane antiferromagnetic spiral in the trans- 
verse field the situation is different (as we discuss later), 
and the structure with rp = tt/2 is realized iZill Below, 
we consider the solution Q12JI. which is stabilized by the 
easy-axis anisotropy, D < 0, for two important orienta- 
tions of the magnetic field. 



A. Magnetic field perpendicular to z-axis 

This field orientation gives rise to the most general 
distortion of the exchange spiral. Deviations of the angles 
satisfying (|13f) are found in the form 



-&sm(2Qr i 



D 



D 



2 ( J o - Jq 



■tg<fo + ¥>cos(2Qr. t ), 



D 



J 3<5 ~ j q J 2Q - Jq 



tg 00, <fi 



D 



J 2Q - J Q 



tg^o- 
(14) 



Note, that according to these expressions, the anisotropy 
causes a small modulation of the spiral winding an- 
gle, (QrJ, which propagates as 2(Qr i ). This natu- 
rally results in appearance of the harmonics at wave- 
vectors ±2Q and ±3Q in the lattice Fourier transform 
of the ground-state spin arrangement, in agreement with 
the early results of Nagamiya et alfi and Cooper and 
Elliott. 15 Another point to be mentioned here is that 
small variations (|14fl . which are treated here in the lin- 
ear approximation, actually contain terms which increase 
with the magnetic field (the most dramatic is the sec- 
ond term in proportional to H 2 ). Therefore, although 
these expressions are valid up to very high fields, they 
fail near the spin-flip transition, at 



jH s ~ 2(J - J Q ) 1 + 



D 



2 ( J 2Q - j q) 



In fact, this failure is indicative of a phase transition 
to a new spin arrangement of the "fan" type, that was 
considered in detail by Nagamiya et aim A simple way 
to study this transition and the behavior of the system 
in the vicinity of the spin-flip field H s , is to expand the 
angles <f>i in small deviations from 7r/2. This approach 
was applied in Ref.ll4l where the field-induced transitions 
in the easy-plane triangular antifcrromagnet CsMnBr 3 
were studied. 

The classical ground-state energy corresponding to the 
spin arrangement of the global minimum found above can 
be easily calculated from Eqs. 0-(lO), 



E GS = NS 2 [Jq + \d-{J -Jq + \d) sin 2 O 



NS J, 



Q 



D 



Nx_l(jH) 2 , (15) 



where we have introduced the transverse spin suscepti- 
bility per site, as follows, 



1 



2(J -Jq)-D 




2(Jo - Jq) \ 2(J - Jq 



0{D 2 ). (16) 



Spin structures are usually studied in neutron scattering 
experiments, by measuring the intensities of the magnetic 
Bragg peaks. At each fe, these are proportional to the 
structure factor, ie the square of the absolute value of 
the corresponding Fourier transform of spin density, S^, 
which, therefore, can be called a Bragg amplitude. Thus, 
for comparison with the experiment, it is important to 
calculate the nonzero Bragg amplitudes which result 
from the solution i|14|) . Introducing the unit vectors e x 



6 



and e y along the x and y axes, respectively, we obtain, 



Sq = e y Ssm(j) 1 + 



D 



2(J - Jq) 



Sq = -(e z - ie x )S cos</> | 1 
-(e 2 + ie x )S cos^o 



= xdlH) 
D 



2(J - Jq 
D D 



S 2 q = e y Ssm(j)o 



^( j 3Q- j q) 2 ( j 2Q - j q 
d 



S 3Q = ( e * _ ie x )S COS 00 



2 ( J 2Q - j q) 

D 



4( j 3 q - Jq) 



-tg 2 </)o 



(17) 



In calculating the Bragg intensities, If, ~ l^fcl 2 , one has 
to keep in mind that vectors Sq and S 3 q have a nonzero 
imaginary part. 



B. Magnetic field along z-axis; low fields 

This case seems to be less complicated, since the dis- 
tortion of the exchange spiral occurs only at low fields, 
H < H c ~ y/\D\ J, and is such that it leaves the spins 
co-planar. This means that <pi — 0, and that the second 
equation in Hll|) and in (| 1 31) is automatically satisfied. 
The first equation is reduced to a rather simple relation, 
from which, to the leading order in D and H, the small 
deviations 69 i are, 







i = ■dsin(2Qr l ) + ?]sin(Qr i ), 

D jH 



V 



J3Q ~ Jq j o + J 2Q - 2J q 



(18) 



It is clear from these expressions, that a leading order is 
insufficient for our approximation. To be consistent to 
the first order in the anisotropy constant D, one has to 
perform the expansions to the second order in 77 ~ H < 
\/\D\ J. Such improved procedure leads to the following 
corrected expression for 



D 



J3Q - j q 



J2Q - j q 1 
J3Q-JQ 4 



(19) 



The corresponding ground-state energy and the parallel 
susceptibility per spin are 



1 



E G s= \S- ( .l n -D- 



Q ' 2 ~ 4( Jq + J 2 Q ~ 2Jq) J ' 
1 



2(J + J 2 q - 2Jq) 



(20) 



so that we can rewrite, r\ = — 2x\\( ; yH). The Bragg am- 
plitudes are also easy to calculate, 



So = e z S 



7# 



X\\jH, 



S Q =7 > ( ez ~ iex ) S 

1 o 



2(J + J 2 q - 2Jq ) 

d + v 2 (J 2 q-Jq) 



1 - 



2 ( J 3Q ~ Jq) 

d + v 2 (J2Q-Jq) 
j sq - j q 



S 2Q = - (e, - ie x )S , 7pr 

4(J + J 2 Q - 2J Q 



- X \\jH(e z - ie x 



S 3Q =( e ^ _ ie x)S 



d + ii 2 (J 2 q-Jq) 
4 ( J 3Q - Jq) 



(21) 



Actually, we can calculate the homogeneous component 
of the spin density, Sq, (it the magnetization) with better 
accuracy than in Eq. (JSTJ, to ~ |D| 3 / 2 . To be completely 
consistent in this order, one should add to 59 i terms like 
/3sin(3Qr,j), with /3 ~ |D| 3 ' 2 . However, retained in the 
first order, these terms would not contribute to So- For 
our purpose it is therefore sufficient to take S9i in the 
form l|18(l . and simply retain more terms in the subse- 
quent expansions. In this way we obtain the following 
correction to So, 



Sq =x\\lH 



1-X\\D 3-4 



+ (mrx\\(j 2Q - J Q ) 



/2Q- 


V 


J 3Q ~ 




(J 2 Q 


- J Q 




- J Q 



(22) 



This expression includes the leading nonlinearity of the 
magnetization, which is often important for comparison 
with experiment. 



C. Magnetic field along z-axis; high fields 

This is the simplest and the best-known case. At some 
critical field, H c ~ y/\D\ J, a usual spin-flop transition 
takes place, and structure becomes the same as that of 
the magnetized exchange spiral, (JSJ), 

Sk = e z S sin ^o^fe^ + + i e y)S cos (^qS^ q 

+ (e x — ie y )S cos 4>o$k -Q • (^3) 
In this case, 



E GS =NS 2 J Q 



4( Jo -Jq+D) 



sin <po 



2( Jo -Jq+DY 



(24) 
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These expressions hold up to a full saturation. The corre- 
sponding spin- flip held, H s , is defined from the condition 
sin^o = 1. Comparing the expressions (1201) and l|24(l for 
the ground state energy, one obtains the following uni- 
versal formula for the spin- flop field H c , 



\ 



2\D\(Jo-J Q ) 1 



Jo - Jq 
J 2Q - J Q 



(25) 



V. EASY-PLANE ANISOTROPY (D > 0) 

For the easy-plane Heisenberg helimagnet, spin reori- 
entation in magnetic field reveals no new features in ad- 
dition to those described in the previous section. All 
spin configurations occurring in this case have their ex- 
act analogs for the easy-axis helimagnet. Therefore, it 
is sufficient to establish the correspondence between the 
cases. The simplest case is that of magnetic field di- 
rected along the z axis. The spin arrangement in this 
case appears to be exactly the same as in Sec. 4.3, and is 
described by the same expressions, Eqs. l (2*5j l and (J52J. 
However, here we must keep in mind that constant D is 
now positive. 

If the field is applied in the easy plane, eg along the x 
axis, the spin arrangement in low fields is similar to that 



described in Sec. 4.2. This correspondence, however, is 
not as direct as the one above. Since all spins lie in 
the same plane, the solution of Eqs. (11) is obtained 
with 6i = 7r/2 and fa = Q r i + 8 fa. The first equation is 
then automatically satisfied, while the second one, for fa, 
becomes exactly the same as equation for 6i in the case 
of a small field along the z-axis in the easy-axis case, 
D < 0, H || z, H < H c , but with D = 0. Since now the 
symmetry breaking in the (easy) spin plane is different, 
we must consider two possible solutions. 

A. First, we assume 8 fa = r)smQr i + r &sm.2Qr i , which 
is exactly the same as considered in Sec. 4.2. In this case 
we find that r\ and d are given by Eqs. (18) and (19), but 
with D — 0. It is clear from the transformation Eq. (7) 
that the resultant Fourier components of the spin density 
are given by the same expressions (21), only with e x and 
e z replaced by e y and e x , respectively. Evidently, such 
solution can be stabilized by the magnetic field in the 
case of a nearly ferromagnetic helimagnet, where spins 
tend to align along the field direction. 

B. Now consider the solution with 8 fa = ir/2 + 
77 cos Qr i — i?sin2Qr i . Parameters 77 and -d are again 
given by Eqs. (18) and (19). As discussed above, in our 
approximation this solution is degenerate with the previ- 
ous (A), and differs from it only by the phase multipliers 
in Bragg amplitudes, 



Sn = 



'2Q 



X\\lH 
1 



l + (^H) 2 x^J 2Q -J Q ) 



J 2Q - J Q 1 



(e y + ie x )S 



1 — 



2 J 2Q-JQ 

2 ( J 3Q - J Q) 



-e v Sr] 2 



1 J 2Q - j q 

2 J sQ - j q 



S 3 q = -(e y + ie x )Sr] 



A{J + J 2 q~2Jq) 2 

J 2Q j q 



-X\\jH(e x -ie y ), 



*( j 3 q- j q)' 



(26) 



It was shown by Chubukov 12 for the case of a six- 
sublattice commensurate helimagnet CsMnBr 3 , that it 
is a type-B structure that is stabilized by magnetic field 
in the sixth order, ~ H 6 . A non-zero frequency, cor- 
responding to the magnon gap ~ H 3 , appears in the 
magnetic resonance spectrum. Magnetic field favors this 
solution for an antiferromagnetic spiral because in this 
case it tends to align all spins perpendicular to its direc- 
tion. In a general case of n-sublattice antiferromagnetic 
spiral, structure is stabilized by the magnetic field 
in the n-th order, as we have discussed in Sec. 4. 

Above the spin-flop transition which occurs for H _L z 
and D > at same H c , given by 1)25(1 . the spin structure 



becomes exactly the same as that discussed in Sec. 4.1 
(assuming H || y), and is described by the same expres- 
sions, but with D > 0. 



VI. CONCLUSION 

In this paper we presented a general pertutbative 
treatment of the spin reorientation process for a weakly 
anisotropic spin spiral described by the classical Heisen- 
berg spin Hamiltonian (1) in a magnetic field. The cal- 
culations were performed up to the first order in the 
anisotropy constant D. Previously, each particular case 
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of this general problem was treated individually, typically 
in the framework of the multi-sublattice model. Many ex- 
amples of such approach are found in the literature de- 
voted to calculating the spin reorientation in "triangular" 
antifcrromagnets on hexagonal Iatticeii2iiii4 Apart from 
the loss of generality, in many important non-trivial cases 
such calculations are very difficult and lead to involved 
formulas which depend on a particular setup of the ex- 
change couplings. As a result, in most cases neither the 
spin-waves spectrum, nor even the magnetic resonance 
frequencies, can be calculated in an explicit analytical 
form on the basis of the sublattice model. In contrast, 
the description presented here is quite general, and the 
results are rather simple. They provide a solid basis for 
spin- wave calculations and for the estimates of the spin- 
fluctuation contribution to the ground-state energy which 
can often be very important (even for selecting the cor- 
rect ground state). Although here we calculated the spin 
configurations only for two principal orientations of the 
magnetic field, the derivation can be easily generalized 
for any field orientation. To this end, one must simply 
retain all relevant terms in Eqs. (11), (13), and admix 
the sin(Qr,) and cos(Qrj) rotations to the angle varia- 
tions 8<bj, 86 j. Thus, for example, some results reported 
in Ref. [lj can be obtained in a more general way. 

As discussed above, one can treat the magnetic struc- 
tures and the long-wavelength spin dynamics within the 
perturbation approach based on exchange approximation 
in the framework of the phenomenological macroscopic 
Lagrangian theory, as was proposed by Andreev and 
Marchenkoii^ For the systems described by the Hamil- 
tonian (1) the leading relativistic corrections are equiva- 
lent to our treatment of the anisotropy. The macroscopic 



theory has the advantage of being less sensitive to the 
fluctuation effects (the mean-field ground state is chosen 
better, to partially account for the fluctuations). How- 
ever, it also has two substantial disadvantages in com- 
parison with the microscopic description presented here. 
First, it is not aimed at (and is not suitable for) calcu- 
lating a microscopic spin-density distribution, or a spin- 
wave spectrum in the whole Brillouin zone. Therefore, 
the results cannot be directly compared with neutron 
scattering experiments which are the most sensitive and 
powerful tool for studying spin systems. Second, to re- 
produce the higher-order effects, such as nonlinear mag- 
netization given by Eq. Il22[| . macroscopic theory requires 
introducing additional phenomenological constants. It is 
clear, that microscopic theory which is developed here in 
the classical-spin approximation is far from being com- 
plete, and requires additional treatment of quantum cor- 
rections. This point, as well as calculation of the spin- 
wave spectra, are beyond the scope of this paper, and 
will be subjects of our further studies. 

Acknowledgments 

We wish to thank L.-P. Regnault, J. Villain, and 
A. F. Andreev for interesting discussions of the prob- 
lem and for illuminating remarks. We also thank 
A. S. Borovik-Romanov, L. A. Prozorova, A. I. Smirnov, 
and S. S. Sosin for reading the manuscript and for use- 
ful comments. The work of I. A. Z. was supported by 
Grant M3K000 of the International Science Foundation. 
M. E. Zh. acknowledges financial support from the Japan 
Society for the Promotion of Science. 



1 A. Yoshimori, J. Phys. Soc. Jpn. 14, 807 (1959). 

2 J. Villain, J. Phys. Chem. Solids 11, 303 (1959). 

3 T. Nagamiya, in Solid State Physics, edited by F. Seitz, 
D. Turnbull and H. Ehrenreich, Vol. 20, 305 (Academic 
Press, New York, 1967); T. Nagamiya, T. Nagata, and 
Y. Kitano, Progr. Teor. Phys. 27, 1253 (1962). 

4 T. A. Kaplan, Phys. Rev. 124, 329 (1961). 

5 D. H. Lyons and T. A. Kaplan, Phys. Rev. 120, 1580 
(1960). 

6 W. B. Yelon and D. E. Cox, Phys. Rev. B 6, 204 (1972); 
Phys. Rev. B 7, 2024 (1973). 

7 M. Eibshutz, R. C. Sherwood, F. S. L. Hsu, and D. E. Cox, 
AIP Conf. Proc. 17, 864 (1972). 

8 I. A. Zaliznyak, Solid State Commun. 84, 573 (1992). 

9 A. V. Chubukov, J. Phys. C 17, L991 (1984). 

E. Rastelli, L. Reatto, and A. Tassi, J. Phys. C 18, 353 
(1985); 19, 19933 (1986). 

1 T. Ohyama and H. Shiba, J. Phys. Soc. Jpn. 62, 3277 
(1993); 63, 3454 (1994). 

2 A. V. Chubukov, J. Phys. C 21, L441 (1988). 

3 H. Tanaka, S. Teraoka, E. Kakekashi, K. Iio, and K. Na- 
gata, J. Phys. Soc. Jpn. 57, 3979 (1988); H. Tanaka, 



Y. Kaahwa, T. Hasegava, M. Igarashi, S. Teraoka, K. Iio, 
and K. Nagata, J. Phys. Soc. Jpn. 58, 2930 (1989). 

14 S. I. Abarzhi, A. N. Bazhan, L. A. Prozorova, and I. A. Za- 
liznyak, J. Phys.: Condens. Matter 4, 3307 (1992). 

15 B. R. Cooper, R. J. Elliott, S. J. Nettel, and H. Suhl, Phys. 
Rev. 127, 57 (1962); B. R. Cooper and R. J. Elliott, Phys. 
Rev. 131, 1043 (1963). 

16 I. E. Dzyaloshinsky, Zh. Eksp. Teor. Fiz. 47, 992 (1964) 
[Sov. Phys. JETP 20, 665 (1965)]. 

17 Yu. A. Izyumov, Usp. Fiz. Nauk 144, 439 (1984) [Sov. 
Phys. Usp. 27, 845 (1984)]. 

18 X. Zhu and M. B. Walker, Phys. Rev. B 36, 3830 (1987); 
M. L. Plumer, A. Caille, and K. Hood, Phys. Rev. B 39, 
4489 (1989). 

19 A. F. Andreev and V. I. Marchenko, Usp. Fiz. Nauk 130, 
39 (1980) [Sov. Phys. Usp. 23, 21 (1980)]. 

20 I. A. Zaliznyak, L. A. Prozorova, and A. V. Chubukov, J. 
Phys. Condens. Matter 1, 4743 (1989). 

21 S. I. Abarzhi, M. E. Zhitomirsky, O. A. Petrenko, S. V. 
Petrov, and L. A. Prozorova, Zh. Eksp. Teor. Fiz. 104, 
3232 (1993) [JETP 77, 521 (1993)]. 



